Draft version February 2, 2008 

Preprint typeset using I^T^jX style cmulatcapj v. 7/8/03 



GIANT MOLECULAR CLOUDS IN M33 
II - HIGH RESOLUTION OBSERVATIONS 

E. Rosolowsky 1 , G. Engargiola, R. Plambeck and L. Blitz 

Department of Astronomy, University of California, Berkeley, CA 94720 
Draft version February 2, 2008 

ABSTRACT 

We present 12 CO(l — > 0) observations of 45 giant molecular clouds in M33 made with the BIMA array. 
The observations have a linear resolution of 20 pc, sufficient to measure the sizes of most GMCs in 
the sample. We place upper limits on the specific angular momentum of the GMCs and find the 
observed values to be nearly an order of magnitude below the values predicted from simple formation 
mechanisms. The velocity gradients across neighboring, high-mass GMCs appear preferentially aligned 
on scales less than 500 pc. If the clouds are rotating, 40% are counter-rotating with respect to the 
galaxy. GMCs require a braking mechanism if they form from the large scale radial accumulation of 
gas. These observations suggest that molecular clouds form locally out of atomic gas with significant 
braking by magnetic fields to dissipate the angular momentum imparted by galactic shear. The 
observed GMCs share basic properties with those found in the Galaxy such as similar masses, sizes, 
and linewidths as well as a constant surface density of 120 Mq pc~ 2 . The size-linewidth relationship 
follows AV oc r °- 45 ± 02 j consistent with that found in the Galaxy. The cloud virial masses imply that 
the CO-to-H2 conversion factor has a value of 2 x 10 20 H2 cm _2 /(K km s _1 ) and does not change 
significantly over the disk of M33 despite a change of 0.8 dex in the metallicity. 
Subject headings: 



1. INTRODUCTION 

Many theories address molecular cloud formation, but 
few observations can distinguish between them. As a 
result, there is little con sensus on which process dom- 
inates l|Elmegreenl 119931 and references therein). One 
avenue for testing these theories is to study molecular 
clouds in a wide range of environments and to determine 
which aspects of the environment set the cloud prop- 
erties. We have adopted this approach in a study of 
the Local Group galaxy M33. lEngargiola et al.1 (2003, 
Paper I) completed a survey of the entire optical disk 
of M33 in 12 CO(J = 1 -> 0) with sufficient sensitiv- 
ity to detect all molecular clouds more massive than 
1.5 x 10 5 M Q . This survey represents the first flux-limited 
survey of giant molecular clouds (GMCs) in any spiral 
galaxy. While there are complete surveys of molecular 
gas in the Milky Way, blending of emission along the 
line of sight make it impossible to obtain a complete cat- 
alog. The only other flux-limited survey of a molecular 
cloud population in another galaxy was completed by the 
NANTEN group st u dying the Large Magelle nic Cloud 
ijMizuno et al.ll200l(K iWilson & ScovTilel l(l990j) used the 
OVRO interferometer to study M33 GMCs at high res- 
olution (20 pc), but their observations were limited to 
clouds in the central 1.5 kpc of the galaxy that are asso- 
ciated with optical extinction and bright single-dish CO 
detections. 

This paper is a high resolution follow-up to the sur- 
vey of Paper I, targeting fields known to contain bright 
GMCs over a large range of galactic radii. The 13" obser- 
vations presented in Paper I have sufficient linear resolu- 
tion (50 pc) to resolve the emission into molecular clouds, 
but not to measure the cloud sizes. The higher resolution 
observations presented here (20 pc) resolve most of the 
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GMCs, and provide additional information about their 
masses, angular momenta and morphologies. These data 
show that the GMCs in M33 are similar to those found 
in the Milky Way. In addition, the data constrain cloud 
formation theories, and they suggest that the CO-to-H2 
conversion factor is constant across the galaxy. 

2. OBSERVATIONS 
2.1. Interferometer Observations 

We observed 17 fields in 12 CO (J = 1 -»• 0) using 
the BIMA array ijWelch et alJ TlQQfi) in Fall 2000 and 
Spring 2001. We chose fields that contain the highest 
mass GMCs cataloged in Paper I (Figure QJ. The obser- 
vations were made in the C-array configuration yielding 
a synthesized beam FWHM of ~ 6". We observed the 
CO line in the upper sideband with a velocity resolution 
of 1.016 km s" 1 . System temperatures typically ranged 
between 420 K and 730 K. Table □ lists the field center, 
synthesized beam size, integration time, and rms noise 
for each field. We tried to obtain < 0.2 Jy beam -1 rms 
noise in each field while maintaining good UV coverage. 
Three additional fields (Fields 18 through 20 in TableQJ 
were observed in conjunction with the survey presented 
in Paper I and are discussed in more detail in that pa- 
per. Because of the longer integration times and smaller 
synthesized beam sizes, the high resolution observations 
can detect lower mass clouds than the survey in Paper I. 

Our data reduction followed the algorith m developed 
for th e BIMA Survey of Nearby Galaxies ijRegan et alJ 
1200 1|) . Only the C-array data were used in producing 
the final data cubes because of their significantly lower 
noise levels. Observations of Uranus (Fall 2000) and 
Mars (Spring 2001) were used to establish the flux of 
the phase calibrator 0136+478 (3.6 Jy in Fall 2000, 2.7 
Jy in Spring 2001). Bootstrapping the quasar flux to the 
M33 data established the flux calibration. Binning veloc- 
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Fig. 1. — The locations of the 20 fields observed by the BIMA 
interferometer and the 18 locations observed by the UASO 12- 
m overlaid on smoothed / band DSS image of the galaxy. The 
locations of the catalog sources from Paper I are plotted as grey 
circles and were used to target the C-array and 12-m observations. 
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Observations of 20 fields in M33 using the BIMA 
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a For a 2 km s 1 channel. 

b These observations of these fields are discussed in Paper I. 



ity channels set the final velocity resolution at 2.00 km 
s _1 . We used uniform weighting of the UV data in the 
inversion, choosing a grid spacing that resulted in 1.5" 
pixels for the final maps, and we cleaned the cube using a 
Hogbom algorithm after correcting for the primary beam 
attenuation. 



2.2. Single Dish Observations 

Since interferometers filter out zero-spacing flux from 
the observed object, complete flux recovery requires ad- 
ditional observations using a single dish telescope. To 
measure the total flux from the GMC fields, we observed 
18 fields over the course of three nights using the Uni- 
versity of Arizona Steward Observatory 12-m telescope 
at Kitt Peak. The selected fields for the 12-m observa- 
tions are displayed in Figure ^ To check for faint emis- 
sion, each field was observed to attain a signal-to-noise 
ratio at the peak of the CO line of at least 10. Every 
two hours, Jupiter or Saturn was observed to optimize 
the pointing and focus of the telescope. During sunset, 
the pointing and focus checks were done before observ- 
ing each new field to compensate for rapid variations in 
the thermal stresses on the telescope. Both polarizations 
were observed with the 500-kHz and 1-MHz filter banks 
in parallel mode. The median system temperature was 
350 K for the observations and the median column of 
water was 0.7 mm. 

The observations used absolute off positions which tar- 
geted low-significance detections in the D-array survey 
(see Paper I). These observations help to establish the 
completeness limit in the survey by determining the frac- 
tion of real detections of low significance. The positions 
of the off fields were selected to be well separated in ve- 
locity from the strong molecular cloud emission in the 
on position. Emission in the off fields appears as absorp- 
tion features in the resulting spectra. Only two of the off 
fields contained emission at the expected velocities. 

We reduced the data by flagging bad channels in ev- 
ery spectrum, combining the polarizations to reduce the 
noise, and fitting a quadratic baseline. The vane calibra- 
tion of the antenna temperature was found to be accu- 
rate to better than 10% using continuum observations of 
Jupiter and line observations of DR21(OH). All spectra 
were observed to an rms noise level of at most 10 mK 
in a 2.6 km s" 1 channel, equivalent to a molecular gas 
surface density of of 0.1 M pc~ 2 . 

Finally, we observed a 5 x 5 pointing mosaic centered 
on the most massive cloud complex in the galaxy (Field 

17 in Tabled]) to map all flux associated with this cloud. 
The single dish map was combined with both the C and D 
array BIMA data using MIRIAD's imcomb routine. The 
difference between the fully-recovered data cube and the 
cube generated solely from the interferometer data shows 
the location and spatial structure of the emission that is 
detected only by the 12-m telescope. 

2.3. Flux Recovery 

By comparing the interferometer and single dish data, 
we estimate the amount of flux that has been resolved 
out by the interferometer or n ot selected by our cloud 
identification techniques f ^2.4[) . We compare each of the 

18 single dish spectra to a representative spectrum gener- 
ated from the interferometer data by averaging together 
the spectra in the data cube and correcting for the beam 
profile of the 12-m. To improve the signal to noise in 
the interferometer spectrum, only pixels with significant 
(> 2er) CO emssion are included in the average. Both 
the C-array observations and the D-array observations 
discussed in Paper I recover less than 100% of the flux as 
expected. The median ratios of the interferometer and 
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single dish integrated intensities {Leo) are: 

(£cO,c/£cO,12-m) = 0.53l°;o9 and 

(Lco,d/ L C OA2-m) = 0-62l ; 13 . 

Because the distribution of this ratio is asymmetric about 
its median, we quote the positive and negative disper- 
sions separately. The 12-m observations indicate that a 
significant fraction of emission is resolved out or masked 
out (see S12.4JI . However, comparing the peak antenna 
temperatures in the spectra shows that: 

(T AfC /T J i,i2_ m >=0.77ig;?? and 

(T Aj n/T A , 12 - m ) =0.89±g:l|, 
while the widths of the interferometer lines are narrower: 
(Au C o,c/Ai;co,i2-m)=0.7li^ and 
(Auoo J D/A«co J i2-m> = 0.7218$. 

Three effects can lead to a discrepancy between the in- 
terferometer and single dish data: (1) calibration errors, 
(2) extended CO emission resolved out by the interfer- 
ometer, and (3) limited sensitivity or masking effects in 
the interferometer data. 

Calibration discrepancies would appear as a rela- 
tive scaling of the observed spectra by a constant fac- 
tor. However, the ratios of the peak temperatures and 
linewidths indicate better flux recovery in the middle of 
the line. We conclude that scaling alone cannot not be 
the source of flux loss. 

Observations of the off positions using the 12-m tele- 
scope indicate that there is no extended CO emission 
on kiloparsec scales with sufficient surface brightness to 
account for the discrepancy between the interferometer 
and the single dish observations. Each of the 18 single 
dish spectra have a 3er noise level of 0.07 K km s _1 (0.3 
Mq pc -2 ); however, accounting for the observed differ- 
ence requires a surface brightness of 1.0 K km s _1 (4.3 
Mq pc~ 2 ). To search for even fainter emission, we av- 
eraged together the 16 empty off positions, shifting the 
velocities to align any faint emission at the H I veloc- 
ity. No line is present in the resulting spectrum, which 
has a 3cr surface brightness sensitivity of 0.019 K km s _1 
(0.08 Mq pc -2 ). Since the off positions are separated 
by projected distances of roughly a kiloparsec, we con- 
clude that a CO emission component on this scale can 
account for no more than 2% of the discrepancy between 
the interferometer and the single dish fluxes. 

The extra emission detected in the single-dish spectra 
is likely from low mass clouds surrounding each GMC ly- 
ing below our clipping threshold. Figure compares the 
map of Field 17 generated from interferometer data alone 
with the map generated by merging interferometer data 
with the single dish map. The greyscale image is the fully 
recovered map and the contours represent the additional 
emission seen only in the fully recovered map. The rela- 
tive levels of the contours to the greyscale are indicated 
at the bottom of the scale bar, showing the additional 
flux recovered is a very low surface brightness compo- 
nent. Taking the Fourier transform of the added emission 
shows that > 80% of the added flux would be detected 
by the interferometer with more integration time. The 
spatial filtering of the interferometer occurs only at UV 
distances of < 5kA and does not account for the dis- 
crepancy between the two telescopes. Most of the added 
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Fig. 2. — Map of emission in Field 17 of the interferometer 
observations. The greyscale shows the map generated by combining 
the interferometer and single-spacing data. The contours show the 
distribution of flux present in the combined map but not in the 
interferometer map. The contours are spaced at 0.5 K km s _1 
starting at 0.5 K km s — approximately la for this map. The 
relative levels of the contours are indicated as lines on the scale 
bar, showing the diffuse emission is of low surface brightness. The 
three diamonds overplotted are the locations of additional clouds 
found in deep C-array integrations. 



emission is excluded from the interferometer map simply 
by virtue of being below the significance threshold for our 
masking methods. Also plotted in Figure are the loca- 
tions of low mass clouds seen only in the C-array data, 
where the synthesized beam is well matched to the cloud 
sizes. Two of these clouds appear in the distribution of 
diffuse emission suggesting clumpy structure within this 
envelope. The diffuse emission seems principally asso- 
ciated with the GMCs, and a discussion of what this 
observation implies for cloud formation appears in t)4.4l 
We conclude that the interferometer observations detect 
the central clouds of complexes in the galaxy containing 
high mass GMCs (M > 1.5 x 1O 5 M , Paper I). 

2.4. Cloud Identification 

We defined clouds in the position-velocity interferom- 
eter data cubes as connected regions in three dimensions 
containing more than 50 pixels, each with a value greater 
than 2a rms . Our criteria guarantee that at least five sta- 
tistically independent regions are included in each cloud 
because interferometer data cubes are oversampled by a 
factor of 3 in each spatial direction, though the veloc- 
ity channels are nearly independent. For five joint 2a 
detections, the likelihood of false detection is 5 x 10~ 4 
in these data cubes which have ~ 10 5 independent ele- 
ments. Analysis of the noise in the data indicates that 
the distribution is sufficiently Gaussian to justify such 
a treatment. We simulated 10 3 UV data sets mimick- 
ing the observation conditions, correlator configuration, 
and baselines of the array. The cloud extraction program 
found no false detections for these selection criteria. 

This algorithm extracts 45 molecular clouds from the 
data cubes, including all clouds in these fields cataloged 
in Paper I. To account for low amplitude signal at the 
edges of the clouds, the regions were expanded by 1 pixel 
in both the position and the velocity directions. Figure 
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shows maps of the individual clouds. Five of the 45 
clouds are not completely recovered since they lie near 
the half power point of the primary beam. Because their 
peaks are included in the maps, we estimate that the 
truncation will not increase their masses and areas by 
more than a factor of 2. 

3. GMC PROPERTIES 

The observed parameters of the 45 clouds are listed in 
Table El The ana lysis uses a distance to M33 of 850 kpc 
l|KennifflitH l i99Sh . 

3.1. Cloud Masses 

We measure the GMC masses using two methods: (1) 
assuming a constant CO-to-H2 conversion factor and (2) 
applying t he virial theorem. T he CO-to-H2 conversion 
factor (X. iBloemen et aTl [T986) converts between inte- 
grated CO line intensity (Leo) and molecular hydrogen 
column density (iV(H2)). We use a conversion factor of 

I.2x 10 2 ° *(H 2 )/cm-» 
Lqo /K km s 

ijStrong k. Mattoxlll996l Dame et alJl2001|) . In convert- 
ing the column density to a cloud mass, we include the 
mass of helium with a number fraction of 9% relative to 
hydrogen nuclei. This gives the luminous mass of the 
cloud: 

Mco .„{ X \ ( L co _^ (A£ 



4.3 



(1) 



M Q ~ \2 x 10 2 V VK km s"V \pc^ 
where A c \ is the projected surface area of the molec- 
ular cloud. Our observations typically have a rm s = 
1.0 K km s" 1 (4.3 M pc~ 2 ). Errors m measuring 
masses using the X factor arise not only from the un- 
certainty in the X factor but also from the flux calibra- 
tion and the noise in the data. Determining the true 
masses of entire GMCs is further complicated by the 
presence of low surface brightness emission around the 
central clouds. 

For the 36 resolved clouds, we can also compute the 
mass using the virial theorem. The virial mass requires 
a measurement of the velocity width of the GMC and the 
cloud size. We used the intensity- weighted second mo- 
ments of the position and the velocity to determine the 
size an d linewidth respectively, following Sol omon et alJ 
l)1987l hereafter S87). The radius of the cloud r is de- 
fined as 



1.36 



2^L 



(2) 



where a a and as are the dispersions in right ascension 
and declination converted to projected distance and <7f, m 
is the beam width. The coefficient scales the result to 
match the definitions of cloud radius given in S87 to fa- 
cilitate comparison with Milky Way GMCs The 
scaling is chosen so that the radius of a circular cloud 
would be the same using either measurement method. 
Our size measure is independent of the orientation of the 
coordinate system in which the moments are calculated. 
With the noise level of our data, this radius measurement 
works best for cloud radii larger than 10 pc, the projected 
beam radius. With this definition of r, the virial mass of 
the cloud is then: 



M- 



VT 



1427 



( °l 



\km s 



pc 



(3) 



This virial mass estimate implicitly assumes a density 
profile of p oc r " 1 . These methods repro duce the virial 
masses found bv lWilson fc Scovillei (HoOflh to within 40% 
for the four clouds that are well-resolved in both samples. 

The virial mass estimate assumes a spherically sym- 
metric cloud so that the velocity dispersion in the radial 
direction corresponds to the size scale in the plane of the 
sky. Therefore, the most reliable virial mass estimates 
come from clouds that appear round on the sky. We 
define the spatial filling fraction / as 



f = 



4A cl 
ir£ 2 

max 



(4) 



which compares the area actually occupied by the cloud 
to the area that the cloud would occupy if it were circular 
with diameter equal to £ max , the cloud's largest linear 
dimension. We restrict analyses involving virial masses 
to the 29 "round" clouds, defined as having / > 0.5, to 
minimize errors in the virial mass estimates. Clouds with 
small filling fractions appear as elongated structures and 
suggesting that they may be decomposed into multiple, 
dynamically distinct GMCs (e.g. Clouds 5, 10, 19, or 21 
in Figure [5J. Figure 0] compares the virial and luminous 
(X factor) masses for the 36 resolved clouds in our study. 
A linear fit between the virial and CO masses for the 29 
round clouds finds M VT =Jl.Q±0.1)M C o - (0.1 ±0.2) x 
10 4 Mq with a reduced x 2 value of 0.9. The slope of 
unity implies that our adopted value for the X factor is 
consistent with the virial masses. 

In >I2.3I we found that our observations underestimated 
the integrated intensity and linewidth of the CO emis- 
sion by 50% and 25% respectively. Recovering more of 
the flux with higher sensitivity observations would in- 
crease the luminous mass as well as the size and linewidth 
used in the virial mass. These increases tend to offset 
each other so the ratio of the luminous mass to the virial 
mass remains roughly constant. To investigate how omit- 
ting 50% of the flux in a field affects our mass estimates, 
we modeled the clouds as azimuthally symmetric bright- 
ness distributions with a power law surface brightness 
(Lco( r ) oc r"). The linewidth scales with the size follow- 
ing the relation found in H4.2I Using the above methods, 
we calculated Mco /Mvt f° r the central 50% of the emis- 
sion and for the entire cloud. For all values of u> between 
—2 and 0, the ratio for the central 50% is never more 
than 1.5 times higher than the ratio for the entire cloud. 
This estimate is an upper limit since some of the emis- 
sion is likely from dynamically distinct, low mass clouds 
surrounding each GMC, implying the interferometer re- 
covers significantly more than half the emission from the 
central cloud. We conclude that there is a good corre- 
spondence between virial and luminous masses derived 
from the interferometer data. 

The values of masses, radii and linewidths, are com- 
parable to those found in the Milky Way (S87). The 
empirical relations between these quantities also match 
those found for the Milky Way molecular clouds and this 
similarity is discussed in detail in H4.2I 

3.2. Velocity Gradients 

The high resolution CO observations allow measure- 
ment of the velocity gradient across the molecular clouds. 
Many mapping studies of Milky Way molecular clouds 
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Fig. 3. — Maps of individual molecular clouds in M33. The greyscale map is the CO brightness, integrated over all velocities contributing 
to the cloud. Each box is 300 pc square. The maps are annotated in the corner with "7V:Lco" where N is the number of the cloud in 
Table 1^1 and Lqq is the peak integrated intensity (K km s _1 ) in the cloud which sets the top of the greyscale. The bottom of the greyscale 
is set at K km s" 1 . Other clouds in the same field of view with higher peak intensities appear saturated. The typical beam is plotted in 
the first panel of the figure. Plotted over the greyscale are contours of integrated intensity for the cloud under consideration. The spacing 
between the contours is 2 K km s" 1 (~ 2cr rms ) if the peak integrated intensity (Lqq) is less than 20 K km s _1 and 4 K km s — 1 otherwise. 
The maps show a wide variety of morphologies and brightnesses. 



note such gradients and are summarized in Phillips 
(1999). The clouds in M33 also show significant velocity 
gradients, with magnitudes comparable to Milky Way 
clouds. To measure the gradient, we derive the veloc- 
ity centroid at each position in the cloud using the in- 
tensity weighted first moment of the masked spectrum. 
Then, we fit a plane to the resulting velocity centroid sur- 
face. The coefficients of the fit define Vu c , the gradient 
in th is surface and iy , the radial velocity of the cloud, 
after iGoodman et al.1 |l993). The vector quantity Vt^ 
can also be expressed as the gradient magnitude, |Vw c |, 
and the position angle of the gradient 4> cloud- The latter 
is measured east from north to facilitate comparison with 
the orientation of the galaxy. Errors in quantities derived 



from the fit are propagated from the original data using 
the covariance matrix. 

The results of the gradient fits are given in Table [21 
Planes are good fits to the velocity centroid surfaces with 
the median value of reduced x 2 = 1-4* 2 ,' 2 . f° r the mea- 
sured gradients. Fitting other functions (like paraboloids 
or radial collapse profiles) to the centroid surface yield 

significantly higher values of x 2 ■ The magnitudes of the 
gradients are displayed in the histogram in Figure[SJ The 
magnitudes are comparable to the typical 0.1 km s -1 
pc" 1 for the c louds found in the Milky Way (|Blitzill993t 
lPhillipsi n"999). In addition, there is no significant re- 
lation between position angle and gradient magnitude. 
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Figure shows position- velocity cuts along the derived 
gradients for four clouds in the data set. Adjacent to 
each plot is a position-velocity cut perpendicular to the 
derived gradient. This figure demonstrates that planes 
are good fits to the data since the cut along the gradient 
shows a linear trend and the perpendicular cuts appear 
constant. 

Most authors assume that the velocity grad i ents are 
the s ignature of cloud rotation jBlit d Il993t [Phillips 
1999). If the gradients are due to rotation, then the 
magnitude of the velocity gradient measures the angu- 
lar velocity vector ft projected into the plane of the sky 
|Vi> c | = | O | sin i . The position angle of the gradient is 
4- 'cloud = 0q + 90°. The gradient therefore provides much 
information about the angular velocity and rotation axis 
of a cloud. We find (1) the linear gradients observed 
are the signature of solid-body rotation. (2) The rotation 
energy for these clouds is < 5% of the cloud binding en- 



ergy. (3) The position angles of the rotation axes show 
marginal alignment with the axis of the galaxy, as is seen 
in Figured (4) There are many clouds with position an- 
gles separated by more than 90° from that of the galaxy. 
(5) If the clouds are rota ting, r oughly 40% of the clouds 
are retrograde rotators f ^4.1.4fl . (6) For the median gra- 
dient magnitude of 0.05 km s _1 pc _1 , the rotation pe- 
riod of a molecular cloud would be 125 Myr, significantly 
longer than the ass umed molec ular cloud lifetimes of 10 
to 30 Myr (Paper I. lmitllT99l . 

We measure the angular momentum of the molecular 
clouds from the gradient fits and the cloud radii. For a 
solid body rotator with a power law surface mass density, 
the specific angular momentum 2 j = J/M = (3\Vv c \r 2 

2 Through the remainder of the paper, we will refer to the specific 
angular momentum as the angular momentum. Any case where the 
true angular momentum, J = jM . is intended will be specifically 
noted. 
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Fig. |21 — Continued 



where (3 is a constant set by the mass profile and r is the 
radius of the cloud. For a constant surface density cloud, 
f3 = 0.4, ranging; fro m 0.33 to 0.5 for different density pro- 
files (|Phillipsl ll999). This analysis adopts = 0.4 ± 0.1. 
The derived angular momentum is listed as j\ in Table 
[21 The alignment of the rotation axis with respect to the 
observer affects the magnitude of the observed gradient. 
The value j\ = j sin i where j is the true angular momen- 
tum of the cloud and i is the inclination to the line of 
sight. Without knowing the inclination, ji will underes- 
timate the true angular momentum. To avoid assuming 
solid-body rotation, we calculated a second measurement 
of the angular momentum 



32 



= y2T A (xi,yj,v k ) [{xi - x ) 2 + {yj - y G f 



1 1/2 



(vk - v ) sin# M • ^2 T A (xi,yj,v k ) (5) 

where xcbZ/o, and ^o, are the centers of the cloud aver- 
aged over position and velocity space respectively and 
i?i j is the angle between the point (xi,yj) in the cloud 
and the gradient judged from the center of the cloud. 
This value can be calculated even when the cloud is too 
small to be deconvolved accurately. Even with this mea- 
surement, a cloud can be oriented with its angular veloc- 
ity vector along the line of sight thereby obscuring any 
signature of rotation. The effects of inclination must ulti- 
mately be dealt with statistically under the assumption 
that the clouds are not preferentially oriented towards 

t he observer. 

iBurkert fc Bodenheimerl |2000, BB00) show that tur- 
bulent velocity fields can also produce linear gradients in 
the observational domain. The required turbulent power 
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Table 2 

Cloud properties for clouds found in M33 



Cloud 
Number 


(a,(5) a 
("," ) 


Alternative 13 
Names 


v 

LSR. 
(km s" 1 ) 


r W rilVl 
(km s" 1 ) 


r c 
(pc) 


V 1 

(10 4 M ) 


oo 
(10 4 Mr) 


|Vu| c 
(0.01 km s -1 pc -1 ) 


4> v 

o 


pc km s~ 


32 

1 pc km s — 1 


1 


19, —18 


EPRB19,WS1,2,4,5 


-168. 


7.7 


27. 


31. 


26. 


8.38 


178 


19.2 


25.3 


2 


29, —36 


EPRB132,WS3 


— 149. 


5.5 






6.7 


2.29 


-64 




2.14 


3 


— 82, —41 


WS33,34 


— 162. 


6.7 






3.6 


9.58 


— 82 




6.18 


4 


71, 100 


EPRB31,WS29 


—215. 


6.1 


22. 


16. d 


4.2 


10.1 


-44 


15.3 


47.4 


5 


61, 119 


EPRB31,WS29 


-223. 


11. 


48. 


100 d 


14. 


13.6 


-60 


93.9 


103. 


6 


81, 92 


EPRB30,WS30 


—214. 


4.9 






2.9 


11.9 


19 




10.3 


7 


— 9, —129 


WS10 


— 151. 


8.4 


13. 


17. 


9.9 


8.36 


143 


4.00 


13.4 


8 


26, —117 


EPRB75,WS6 


— 158. 


4.9 






2.9 


3.49 


— 33 




2.44 


9 


27, —129 


EPRB75,WS8,9 


-152. 


4.4 


13. 


4.6 


5.3 


2.20 


—9 


1.06 


4.22 


10 


113, 115 


EPRB18,WS27,28 


— 221. 


6 . 5 


33. 


27. d 


12. 


1.90 


— 16 


6.43 


8.87 


11 


119, 69 


EPRB3,WS21,22,23(?) 


—210. 


9 . 6 


32. 


56. 


41. 


9.53 


24 


30.0 


46.0 


12 


52. —113 


EPR.B16 


— 168. 


7.6 






6.7 


10.6 


32 




12.8 


13 


118, 133 




— 208. 


7.2 


22. 


22. 


11. 


3.29 


-46 


5.06 


4.56 


14 


25, -157 


EPRB75,WS7 


-153. 


4.2 


6.7 


2.3 


3.7 


7.70 


— 91 


1.04 


12.6 


15 


-127, -22 


EPRB4.WS32 


-166. 


10. 


31. 


01. 


45. 


9.57 


-67 


27.5 


45.6 


Hi 


— 138, -57 




-165. 


5.5 


13. 


7.3 


7.0 


1.25 


119 


0.617 


1.64 


17 


153, —58 


EPRB46,WS16,17,18 


-182. 


8.2 


3.8 


4.8 


12. 


14.0 


-4 


0.624 


16.2 


IS 


249, 145 




—224. 


6.0 






3.4 


4.70 


— 176 




3.08 


m 


-4, -340 


EPRB20 


— 134. 


4.8 


49. 


22. d 


23. 


0.685 


52 


4.99 


0.0282 


20 


233, —26 


EPRB7 


— 195. 


6.0 


14. 


9.4 


15. 


8.74 


-17 


4.96 


10.5 


21 


—99, —386 


EPRB22 


— 119. 


5.8 


42. 


27. d 


12. 


4.40 


137 


24.2 


20.4 


22 


291, 144 


EPRB58 


-210. 


4.6 


10. 


4.0 


5.5 


6.75 


93 


2.04 


10.4 


23 


290, -28 


EPRB72 


-196. 


5.1 






6.2 


7.77 


-123 




11.2 


24 


211, 490 


EPRB10 


— 258. 


4.2 


14. 


4.7 


6.7 


4.05 


—5 


2.47 


4.75 


25 


201, 496 


EPR.B10 


—256. 


0.7 


18. 


15. 


15. 


7.24 


— 109 


7.26 


13.9 


26 


— 136, 358 


EPRB15 


—213. 


7.2 


7.8 


7.7 


8.5 


11.9 


— 140 


2.24 


12.0 


27 


— 129, 378 


EPRB15 


—216. 


0.0 


16. 


11. 


5.5 


1.54 


112 


1.16 


4.96 


28 


256, —199 


EPRB6 


— 160. 


. 5 


32. 


25. 


27. 


3.49 


—24 


10.7 


6.92 


20 


330, —19 


EPR.B21 


— 191. 


7.4 


38. 


40. 


21. 


8.91 


— 170 


40.0 


47.0 


30 


— 129. 395 


EPRB15 


—220. 


7.3 


2fc 


26. 


16. 


9.44 


— 157 


19.3 


33.9 


31 


103, 547 


EPR.B2 


—243. 


6.5 


29. 


23. 


32. 


3.58 


— 86 


9.06 


17.6 


32 


132, 559 




-244. 


5.2 






5.8 


4.76 


63 




4.31 


33 


268, 526 


EPRB74 


— 246. 


8.5 


28 


38. d 


8.0 


5.35 


— 55 


12.4 


18.7 


34 


211, 562 




— 253. 


7.4 


8.8 


9 . 2 


7.5 


13.8 


69 


3.29 


16.7 


35 


90, 562 




-246. 


5.4 


7.9 


4.3 


6.2 


9.35 


-135 


1.77 


19.1 


36 


115, 579 




-253. 


6.7 


47. 


41. d 


24. 


3.96 


-54 


26.9 


45.8 


37 


238, 570 


EPRB1 


-248. 


9.6 


38. 


67. 


78. 


7.37 


24 


33.3 


27.6 


38 


218, 625 


EPRB108 


-249. 


6.8 


16. 


14. 


7.1 


9.67 


116 


8.03 


14.6 


39 


295, -352 


EPRB5 


-157. 


9.0 


21. 


32. 


32. 


7.47 


-110 


9.61 


29.3 


40 


548,413 




-248. 


6.2 


13. 


9.3 


9.6 


7.93 


-144 


3.85 


12.8 


41 


552, 431 


EPRB8 


-242. 


8.7 


19. 


27. 


28. 


6.97 


175 


7.66 


11.4 


42 


564, 403 


EPRB9 


-222. 


9.3 


23. 


38. 


30. 


8.22 


-130 


13.7 


28.3 


43 


-362, -821 




-74.0 


6.3 






4.6 


9.95 


-104 




7.18 


44 


618,48 


EPRB11 


-202. 


0.4 


18. 


14. 


6.7 


10.8 


-30 


10.7 


23.4 


45 


627,65 


EPRB11 


-204. 


9.2 


23. 


38. 


25. 


8.56 


32 


14.4 


19.2 



a Right ascension and declination are given in seconds of arc relative to the galactic center at 02000 = l h 33 m 50?8 and (52000 = 
30° 39' 36" 7. 

b Names are those found in Wilson & Scovillc 1 1990), denoted WSXX where XX is the cloud number in Table 2 of their paper 
and those found in Paper I, denoted EPRBXX where XX is the number in that paper. Clouds referred to twice are clouds 
resolved into two separate clouds by higher angular resolution. 

c Properties are derived using methods described in text. Approximate errors are <5r = 5 pc, <5My<p =20%, 8Mqq=40%, 
<5|Vt;|=0.02 km s" 1 pc" 1 ., <5</>=10°, <5ji=0.03 pc km s" 1 . 

d Filling fraction (/, Equation[|J is less than 0.5, so the virial mass estimate is unreliable. 



spectrum, P(k) oc k n , has a value of n between —3 and 
—4, so that the largest size-scales have the most power, 
appearing as gradients. Using their scaling for the mag- 
nitude of the velocity gradient as a function of size gives 
gradients on order 0.08 km s^ 1 pc -1 for the median cloud 
radius of 20 pc, comparable to the observed values. If 
the gradients are due to turbulence and not rotation, 
our measurements change only in that the angular mo- 
mentum derived from the data will overestimate the true 
value of the angular momentum by a factor of 2 to 3 



(BBOO). The orientation of the net angular momentum 
vector is still represented by the position angle of the 
gradient. 

4. DISCUSSION 

4.1. Angular Momentum Constraints on CMC 
Formation 

The measured angular momenta are significantly less 
than expected from simple theories of cloud formation. 
In the absence of external forces, the angular momentum 



GMCs in M33 



9 



lCr 



o 



10" 



• Circular Clouds 
o Elongated Clouds 



•a . 



10 b 

X factor Mass 



1CT 



10 4 M 



Fig. 4. — Comparison of virial and luminous masses. Elon- 
gated clouds are those with a spatial filling fraction less than 0.5. 
Errors for M VT are ±30% and for M co are ±40%. The solid 
line indicates the locus where Mvt = ^CO- Clouds with a spa- 
tial filling fraction greater than 0.5 show good agreement between 
their virial masses and the masses derived from the X factor with 
X = 2 X 10 20 H 2 cm" 2 /(K km s" 1 ). 
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Fig. 5. — Histogram of gradient magnitude values for clouds in 
M33. Negative values are given to those clouds that have a position 
angle differing from the galaxy by more than 90°. The gradient 
magnitudes are comparable to typical values found in the Milky 
Way. Moreover, the magnitudes of the gradients are comparable 
among clouds independent of alignment with the galaxy. 



of a molecular cloud should be equal to the angular mo- 
mentum of the gas from which it formed. In a shearing 
galactic disk, the initial angular momentum for a forming 
cloud is (see Appendix lAl : 



AR 2 . 



(6) 



R=R C 



In this equation, R is the distance from the center of 
the galaxy, R c is galactocentric radius where the cloud 
is found, and V(R) is the galactic rotation curve (from 
iCorbelli &: Schneider! Il997). The radial accumulation 
length AR is half the radial extent of the parent atomic 



Fig. 6. — Examples of four fitted gradient for molecular clouds 
observed with C-array. The left-hand column contains plots of the 
velocity centroids denoted by lcr error-bars from a position-velocity 
cut along the direction of the gradient. The dotted line is a linear 
fit to the gradient. The right hand column is identical to the left 
column, except the position-velocity cut is made perpendicular to 
the gradient. The minimal curvature in the left column and zero 
slope in the right indicate that planes are good fits to the velocity 
centroid surfaces. 



gas. The radial mass distribution of the collapsing gas 
determines the parameter rj. 

In the simplest case, the radial accumulation length 
is set by requiring the mass of atomic hydrogen in a 
cylinder stretching to infinity in the z direction and of 
radius AR to equal the mass of the molecular cloud: 
E H , • TrAi? 2 = M G mc- We adopt 77 = 1/4 in Equa- 
tion E| calculated for the collapse of a cylinder. Using 
the measured masses of molecular clouds and the lo- 
cal value of E HI , we find that in all but 3 of the 36 re- 
solved molecular clouds, the predicted angular momentum 
is higher than the observed value, on average by a fac- 
tor _o£_5^_U^_£ajcidated^he_va2ues of E HI using the maps 
of IDeul fc van der HulstJ l|1987|) a fter recalibrating this 



map t o match the surface density of lCorbelli fc Schneider! 
( 1997). If the velocity gradients are due to turbulence in- 
stead of rotation, BB00 argue that the observed angular 
momentum should be scaled down by a factor of 2 to 3, 
which only increases the discrepancy between the pre- 
dicted angular momentum and the observed value. Sta- 
tistically, projection effects can account for only 20% of 
the discrepancy. 
We refine this simple model as follows: (1) by assum- 
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Fig. 7. — Histogram of the position angles of the GMC rotation 
with respect to that of the galaxy. Clouds show slight preferential 
alignment with the orientation of the galaxy. The line represents 
the predicted distribution for a value of a = 0.69 in Equation 151 
Rough separations between Prograde and Retrograde clouds are 
shown, but the unknown inclinations of the clouds prevent this 
separation from being completely determined. 



ing that the gas far from the galactic plane does not 
contribute to the forming molecular cloud and (2) by in- 
cluding some estimate of projection effects. Instead of 
the progenitor material extending infinitely far from the 
galactic plane, we assume the gas is accumulated over a 
cylinder with its height equal to its diameter, centered 
on the location of the GMC. The three-dimensional dis- 
tribution of atomic hydrogen is unknown, but we assume 
a scale height of H =100 pc throughout the galaxy with 
p(z) oc sech 2 (z/2H). Changing this value by a factor of 
2 does not alter the results significantly. Again, the di- 
mensions of the cylinder are determined by requiring its 
volume to contain a mass of atomic gas equal to that of 
the molecular cloud. Using the radius of the cylinder as 
the radial accumulation length AR, the predicted angu- 
lar momentum is determined from Equation [S] We also 
assume the clouds are solid body rotators with random 
orientations and scale the predicted angular momentum 
down by a factor of (sini) = 7r/4. If, instead, we assume 
the clouds are strictly aligned with the galaxy, the correc- 
tion only differs from ir/4 by 3%. We compare this model 
to the measurements in Figure |H1 The measured angular 
momenta are plotted as ±1<t error bars as a function of 
galactic radius. Figure |H1 shows that the observations are 
always less than the predictions of this Solid Body Model 
(filled circles). The discrepancy is more than an order of 
magnitude for 20 of the clouds. 

4.1.1. The Toomre and Parker Instabilities 

Recent work on star formation at galactic scales sug- 
gests that star formation pre ferentially occurs where th e 
gas disk is Toomre unstable ijMartin fc Kennicuttfe OOl'). 
We calculate the radial accumulation len gth if the clouds 
actually formed from the instability. The lToomrd (1964) 
instability is a gravitational instability in a shearing, thin 
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Fig. 8. — Specific angular momenta (ji) imparted by different 
collapse mechanisms as a function of galactic radius. The measure- 
ments are represented by ±l<r error bars. The two curves are the 
imparted angular momenta as a function of R ga l for the Toomre 
and Parker Instabilities. The filled points are the predicted values 
of the specific angular momentum for the solid body rotator dis- 
cussed in the text. Nearly all measurements of the specific angular 
momentum lie well below the predictions for all models. 



fluid disk that bunches the gas into concentric rings, pro- 
vided the disk is unstable according to the criterion 



nGY,(R) 



< 1 



(7) 



where k is the local epicyclic frequency, o~ v is the (three 
dimensional) velocity dispersion of the fluid, and S(i?) 
is the local surface mass density of the fluid under con- 
sideration. The radial accumulation length is set by the 
original radial extent of the contracting rings: 

ARt _ 2fonm _ v 



where At is the most unstable scale for Toomre collapse. 
The rotation curve determines the epicyclic frequency for 
the galaxy 

'V^_ VdV' 

#2 + RdR i 



k 2 = 2 



We assume that there is sufficient mass in the annulus 
to form molecular clouds which is true over most of the 
galaxy. 

Figure |H1 also shows the predicted values of the cloud 
angular momentum for GMCs formed using the Toomre 
instability. These values only agree with the observed 
angular momenta for a few clouds in the inner 500 pc of 
the galaxy, where the disk is most stable according to the 
criterion in Equation Since all but two of the clouds 
disagree with the predictions of the Toomre instability by 
a large margin, it seems unlikely that the clouds form by 
the action of this instability alone. Moreover, the Toomre 
instability predicts that cloud angular momentum should 
increase with galactocentric radius, which is not seen in 
Figure 03 

The Parker instability is another potential avenue 
of molecular cloud formation llMouschovias et all 119741 
IElmegreenlll982t IHanawa et alJH992j) . The Parker insta- 
bility is a magneto-hydrodynamic (MHD) effect in gas 
coupled to a buoyant magnetic field. As the instability 
evolves, a section of the field bulges out of the plane of 
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the galaxy, causing the atomic gas to flow down the grav- 
itational potential well along the edge of the bulge and 
collect into molecular clouds. The characteristic length 
scale of the instability is Xp « ttH, where H is the scale 
height (|Oiz fcShul 11993ft . The values for the predicted 
angular momentum are shown in Figure |B1 for a 100 pc 
scale height. These values have been corre cted for the 
pitch angle of the magnetic field (6 = 60°, lBecldl2000ft 
since the collapse occurs along the magnetic field lines, 
so that ARp — (Apsin6*)/2. The Parker instability also 
gives a large discrepancy between the observed and pre- 
dicted values. Matching the observed values would re- 
quire a gas scale height of ~ 20 pc, which would imply 
an unrealistically large stellar mass density. 

4.1.2. The Angular Momentum Problem 

We conclude none of the simple theories produce the 
observed values of angular momentum in the molecular 
clouds. Figure IS] compares the theoretical and observed 
angular momentum for the resolved molecular clouds. 
Even using the moment based measurement of angular 
momentum ( j'2 ) , we find an average discrepancy of more 
than a factor of 5 between observation and theory. This 
discrepancy shows that nearly all clouds suffer an angu- 
lar momentum problem independent of assuming solid- 
body rotation. The failure of the Solid Body Model is 
particularly troubling since this model gives the mini- 
mum accumulation length required to gather the mass 
of the GMC from the atomic hydrogen. This constraint 
will limit even the most sophisticated models of cloud 
formation. 

There are at least two possible solutions to the angular 
momentum problem. One is to gather material for molec- 
ular clouds from similar galactocentric radii, thereby 
minimizing the shear in the atomic gas, as may happen 
in th e swing amplification instability ijKim fe Ostrikerl 
12001ft . If clouds have small radial accumulation lengths, 
they must have long azimuthal accumulation lengths. 
Agreement with the observed angular momenta requires 
AR < 30 pc. To form a cloud of 10 5 M Q out of atomic 
hydrogen with a typical surface density of 7 Mq pc -2 re- 
quires azimuthal accumulation lengths of ^500 pc. With 
the maximu m non-circular velocities in t he atomic gas < 
10 km s -1 i)Deul fc van der Huistl 11987ft . this implies a 
limit on the timescale of cloud formation of Tf orm > 50 
Myr, larger than the 10 — 20 Myr cloud lifetimes implied 
in Paper I. This timescale is a significant fraction of the 
time between the spiral arm crossings f or most of th e 
galaxy away from the co-rotation radius l|Newtonl fl980'). 
Consequently, such a long time scale for cloud formation 
has difficulties explaining preferentially finding clouds in 
spiral arms of galaxies. 

A second solution is to brake the material with external 
forces such as magnetic fields. Magnetohydrodynamic ef- 
fects seem to be critical for shedding angular momentum 
in the star formation process and may also act in the 
early stages of cloud formation. If the magnetic field is 
so dynamically significant as to produce an instability, it 
may also provide en ough tension to brake the clouds dur- 
ing their formation. iBeckl l)2000D measures the mean field 
strength in M33 as 6±2/jG, which gives an Alfven speed 
of ^6 km s _1 (assuming a sech 2 density profile to cal- 
culate the midplanc volume mass density). This value 
is comparable to the sound speed in the warm neutral 



medium of most galaxies, suggesting that MHD braking 
and cloud formation via instabilities occur on compa- 
rable timescales. The braking time is set by the time 
needed for an Alfven wave to sweep through a volume 
of diffuse gas with a moment of inertia e qual to that of 
the forming cloud l)Mestel fc Parisl lT984'l. If the atomic 
gas has a constant density, then the Alfven wave must 
pass through a distance comparable to the accumulation 
length for the molecular cloud. Instabilities cause cloud 
formation on a time scale comparable to the sound cross- 
ing time for the accumulation length. Since the crossing 
times for sound and Alfven waves are comparable, we 
conclude that MHD braking is a viable option to slow the 
rotation of forming clouds. Because only 10% of clouds 
show angular momentum approaching the expected val- 
ues from any of the formation mechanisms (e.g. Figure 
121), the MHD braking must occur in the atomic gas or 
during a small fraction of the molecular cloud lifetime. 
We return to the implications of these observations for 
cloud formation in A4.4I 

4.1.3. Angular Momentum Magnitude Correlations 

The molecular clouds in this sample suggest a correla- 
tion between specific angular momentum and mass. Fit- 
ting a power-law relationship between j\ and A/co finds 
j(M) oc M°- 6±01 , though there is significant scatter in 

the data (x 2 = 6.1, Figure |5J. A si milar correlatio n ap- 
pears in the Milky Way data where IPhillipsI <)1999ft finds 
j oc Af ' 7 . If the molecular clouds form by accumula- 
tion of atomic gas, then the angular momentum should 
scale as the accumulation length squared (Equation EJ: 
j oc AR 2 . The mass of the cloud is M = 7rAi? 2 S~^, 
so a correlation is expected with j oc M/E HI . The vari- 
ations in the mean hydrogen surface density are small 
compared to the order of magnitude spanned in mass 
( H4.2JI . so the angular momentum is expected to scale 
with mass for any accumulation model. That j oc M ' 6 
instead of j oc M may be due to the action of the brak- 
ing mechanism, though j oc M is not strongly ruled out 
by the poor quality of the fit. Such correlation between 
j and M is also expected for a turbulent velocity field. 
BB00 find that the velocity gradient scales with clouds 
size r as |Vu| oc r~ 5 so j\ oc |Vi>|r 2 oc r 15 . The clouds 
in our sample have a constant surface density regardless 
of cloud size ( ^4.21) . Therefore, the mass of a cloud is 
proportional to its area M oc r 2 and j oc M 75 , close to 
the observed trend. 

There is no significant variation of specific angular mo- 
mentum as a function of galactic radius, contrary to the 
predictions of all the accumulations mechanisms consid- 
ered here (Figure |HJ- This implies that the angular mo- 
mentum of the clouds does not depend on any of the 
properties that change across the galactic disk, including 
galactic shear, surface mass density of stars and atomic 
gas, interstellar pressure, interstellar radiation field, or 
metallicity. 

4.1.4. Position Angle Correlations 

The position angle of the angular momentum vector 
can be measured from the velocity gradients as discussed 
in ^13.21 The data show correlations between the cloud 
position angles and that of the galaxy as well as among 
the individual GMCs. If the angular momentum is set 
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Fig. 9. — The correlation between the magnitude of the specific 
angular momentum and the mass of the cloud. Only clouds that 
are well resolved with low filling factors are included in this fit. 
The solid line is a least-squares, power-law fit to the data with 
ji oc JV/O 'SiO' 1 accounting for errors in both directions. 



entirely by the galactic shear, then the rotation axes of 
the GMCs should be aligned with the galaxy. The dis- 
tribution of gradient position angles is shown in Figure 
[7| including rough divisions between prograde and ret- 
rograde rotators (assuming the gradients represent rota- 
tion). Because of the unknown inclination of the clouds 
relative to the line of sight, we analyze the distribution of 
position angles statistically. We assume that the clouds 
are distributed with orientations on the sphere set by the 
probability distribution function (PDF): 



1 

2na 



cos ( 



1 



-1 l/o-l 



(9) 



where 9 is the angle between the angular momentum vec- 
tor of the cloud and that of the galaxy, a is a parameter 
that determines the degree of alignment between these 
vectors, and ip is the azimuthal angle. For a = 1, the ori- 
entation angles are randomly distributed on the sphere 
so there are an equal number of prograde and retrograde 
rotators. For a w 0, the clouds are closely aligned with 
the galaxy (a — 0.2 implies that 91% of the clouds are 
prograde rotators). This PDF is motivated solely for 
ease of simulation so that a single parameter a deter- 
mines how well aligned the clouds are with the galaxy. 
By determining the value of a that most closely rep- 
resents the observed distribution of position angles, we 
measure the degree of alignment between the clouds and 
the galaxy. We generated position angle distributions 
following this PDF using Monte Carlo simulations and 
then used a one-sided K-S test to measure the difference 
between the observed and simulated distributions. After 
testing a range of a, we found the most probable value 
is a = 0.69 (with K-S likelihood P = 0.99), implying 
61% of the clouds are prograde rotators. The likelihood 
of a random distribution on the sphere (a — 1) is sig- 
nificantly lower (P = 0.60) implying some alignment be- 
tween the clouds and the galaxy. Nevertheless, the wide 



Cloud Separation [kpc] 

Fig. 10. — The correlation between the position angle of molec- 
ular clouds and their separation for clouds with M > 8 X 1O 4 M0. 
There is a trend that clouds at small distances tend to have their 
velocity gradients aligned with each other. The mean in each bin 
is weighted according to the uncertainty in the measurements and 
the error bars are the error in the mean. The first bin is the average 
of 21 pairs of position angles. 



range of position angles indicates any angular momen- 
tum imparted by galactic shear is randomized in some 
fashion. 

In addition to the cloud-galaxy correlation, nearby 
clouds show a higher degree of correlation between their 
gradient position angles. Figure ^] shows the differ- 
ence in position angle between pairs of clouds as a func- 
tion of cloud separation. This figure uses clouds with 
M > 8 x 10 4 M Q since these clouds have the most signifi- 
cant gradients. If the clouds are randomly oriented with 
respect to each other, the value of this statistic should 
be 90°, but clouds with separations smaller than 500 pc 
show a significantly smaller value implying their angular 
momentum vectors are aligned. This alignment would 
be expected if these neighboring clouds formed from the 
same parent cloud of H I. If clouds are the product of 
large scale motions in the diffuse ISM, this 500 pc length 
may represent the scale of the flows that give rise to the 
clouds. 

4.2. Larson's Laws in M33 

The clouds in the M33 observations have similar radii 
(~ 20 -> 50 pc), masses (~ 10 4 -> 10 6 M ) and velocity 
dispersion FWHMs (5 — > 10 km s _1 ) as those found in 
the Milky Way by S87. We tested these similarities in 
more detail and found the GMCs in M33 were indistin- 
guishable from those in the Milky Way. We emphasize 
comparison of our results with the work done by S87 
because both studies use equivalent methods for deter- 
mining the properties of the molecular clouds and the 
sensitivities are comparable: ~0.7 K km s -1 for S87 vs. 
0.6 — > 1.3 K km s _1 for this work. We only use the 
23 clouds that are well-resolved (r > 10 pc) and round 
(/ > 0.5) in these comparisons since the derived proper- 
ties of these clouds are the most reliable. 

The results of Paper I show that the mass spectrum 
for GMCs in M33 differs from that found in the Milky 
Way: dN/dM oc M~ 2 - 6 for M33 vs. dN/dM oc M -1 - 5 
for the Milky Way (S87). Despite differences in the mass 
distribution, clouds in M33 show the same power law re- 
lationships between size, mass and linewidth as are seen 
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in the Milky Way (Larson's Laws, after iLarsonl |l981), 
see also S87). We plot the size-linewidth and the mass- 
linewidth relationship for the clouds in M 33 in Fig ure ITT] 
along with the data fro m the Milky Way ijSolomon et al.l 
ll987HHever et al.ll2001fl . Where feasible, the cloud prop- 
erties of the Milky Way data have been recomputed using 
the methods in Sj2|to ensure a common basis for compar- 
ison. Fitting power-laws to the two relationships shows 
AV oc r - 45 ^ 02 and M oc AV 4 - 2±0 - 3 . The fits are con- 
sistent with the relations seen in the Milky Way (S87) in 
both normalization and scaling. However, cloud sizes in 
M33 are systematically smaller by ~ 30% in M33 com- 
pared to S87, which is most likely due to differences in 
the size measurement technique. 

Figure IT2l shows the mean surface mass density within 
molecular clouds as a function of cloud mass. This figure 
shows no significant scaling of surface density with mass. 
Similar plots of surface density as a function of cloud ra- 
dius and galactic radius also show no scaling. The slopes 
of linear fits between the mass density and cloud mass or 
cloud radius are consistent with zero. In addition, two 
statistical tests of correlation, the Spearman rank order 
and the Kendall r tests show no significant correlation 
betw een mass, radius and surface density fsee lPress et alJ 
1992). We conclude that GMCs have a surface mass den- 
sity of 120 M Q pc" 2 (2iV(H 2 ) = 6 x 10 21 cm~ 2 ) with a 
dispersion of 60%, which does not vary with cloud mass 
or radius. This constant mass density is also seen for the 
GMCs in the Milky Way HS'.h ll993l). 

The consistency of Larson's laws between the Milky 
Way and M33 shows that the macroscopic properties of 
a GMC (r, AV, N(H.2)) are the same in both galaxies for 
a given cloud mass. Since these macroscopic properties 
follow similar trends, it is reasonable to expect that the 
internal properties of GMCs depend only on the mass of 
the cloud. The GMCs through out the Local Group also 
follow Lar son's Laws; the LM C l)Mizuno et al.l200ifi . An- 
dro meda llSheth et al J 1200 01. and dwarf elliptical galax- 
ies llYounJ2001l 12002) all show some evidence that cloud 
properties are set by cloud mass. Some unifying mech- 
anism must establish the macroscopic properties of the 
cloud solely in terms of the cloud mass. Since turbulence 
could provide a relationship between size and linewidth, 
the balance between turbulent support and self-gravity 
would naturally set the macro scopic properies of a cloud 
solely in terms of its mass l)Larsonl 119811 lElmegreenl 
1989). 



4.3. The CO-to-H2 Conversion Factor 

In M33, we are presented with a unique opportunity 
to study the effects of galactic environment on the CO- 
to-H2 conversion factor (X) since there are no system- 
atic differences due to inhomogeneous observation tech- 
niques or multiple galaxies. Many studies report a vari- 
ation in the X facto r as a function of metallicity e.g. 
lArimoto et alJ (|199(tD and references therein. There is no 
evidence for such a trend in M33. In Figure we plot 
the ratio of the luminous to virial mass Mco /AfyT as a 
function of metallicity for the 23 well-resolved (r > 10 
pc), round (/ > 0.5) GMCs in our study. The metal- 
licity at each point is obtained using the position of 
the cloud and the metallicity gradient (d[0/H]/dR) of 



IHenrv fc Ho^aTdl Ip9l . Fitting a line to the data gives: 
Mm 

— — = (0.03 ± 0.23) ([O/H] - [O/H] ) + (0.88 ± 0.09) 

MyT 

accounting for errors in both the metallicity (0.1 dex) 

and the mass ratio ijPress et all 1992'). The value of \ 2 

for the fit is 0.9. Th e value of x 2 using the relation- 
ship of lArimoto et al.1 1)199(1 dotted line) as a fit gives 
3.31, showing poor agreement (P = 10~ 5 ) with their 
quoted trend. The slope of our fit is consistent with no 
variation in the conversion factor over a range of 0.8 
dex in metallicity (a factor of 6). iBrand fc Wouterlootl 
( 1995) note that the ratio between virial and CO masses 
does not change systematically with radius in the Milky 
Way. In A3. II we argue that interferometer observations 
overestimate the ratio Mqo /A/vt by a t most a factor of 
1.5. This cannot acc ount for the discrepancy with the 
lArimoto et al.1 l)1996|) data, which we attribute to the 
difficulties in synthesizing a homogeneous data set from 
multiple observational studies using a variety of analysis 
techniques. 

The studied GMCs range in galactocentric radius 
from 170 pc to 4000 pc, and there are significant 
changes in galactic properties over this range. The 
interstellar radiation field changes value by nearly an 
order of magnitude over this r ange in galactic radius 
ijWalterbos fc Greenawalj E?96) . Other variations in- 
clude changes in the cosmic ray flux and the midplane 
pressure of the gas. Such robustness of the X factor as a 
GMC mass measure is observed in the Milky Way (S87) , 
but variations among galaxies may still produce differ- 
ent values in the conversion factor. For example, the X 
factor may change for GMCs in even lower metallicity 
systems such as the SMC that have an average metallic- 
ity below the minimum 8.1 observed in M33. 

4.4. Implications for Cloud Formation 

This study of M33 suggests several attributes of the 
molecular cloud formation process. (1) Molecular clouds 
form directly from atomic hydrogen rather than the ag- 
glomeration of small molecular clouds. (2) The formation 
process makes both low mass and high mass molecular 
clouds in regions of cloud formation. (3) The formation 
events are local, i.e. length scales less than 500 pc. (4) 
The progenitor gas has its initial angular momentum dis- 
sipated by MHD effects. (5) The similarity of molecular 
clouds across many environments implies that the molec- 
ular mass distribution sets the properties of molecular 
clouds and star formation. 

Our results complement those of Paper I where it is 
argued that molecular clouds form directly from atomic 
hydrogen. In that paper, the striking correspondence 
between molecular clouds and dense atomic filaments 
strongly suggests cloud formation from atomic gas. This 
work strengthens that suggestion since we find no ev- 
idence for diffuse molecular gas that could serve as a 
precursor to GMCs ( H2.3|l . Additionally, the atomic gas 
is significantly more massive than the molecular compo- 
nent on a global as well as a local (~ 100 pc) scale. If 
diffuse molecular gas were the precursor material, the 
angular momentum problem would be greatly exacer- 
bated, scaling predicted angular momenta up by a factor 
of Shi/Sh 2 . diffuse > 20. We also find that the flux not 
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Fig. 11. — The size-linewidth and mass-linewidth radius relationship for the molecular clou ds in M33. In bot h plots, the overlaid 
grey dots represent the molecular clouds from a Milky Way catalog generated by merging S87 and lHever eTaT] <200ll) . The power-law fits 
give Al/ oc r 45±0 02 an( j m x ^y4.2±0.3 ip ne ra( ji us measurements for the Milky Way have been recalculated to match the methods 
discussed in J3] though a slight offset between the data sets persists. 
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Fig. 12. — Plot of molecular surface density averaged over 
GMC as a function of mass for 23 well-resolved, round clouds. The 
surface density is calculated using the deconvolved radius r and 
the luminous mass Mqo- There is no significant scaling with mass. 
Similar plots as a function of cloud radius and galactic radius also 
show no scaling. The mean surface density (dashed line) is 120 Mq 
pc~ 2 with a scatter of 60%. 



Fig. 13. — Variation of Mco/Myt as a function of metallicity 
for 23 well-resolved, round clouds in M33. A linear regression shows 
no significant effect from changing metallicit ies, shown by the solid 
line. The dashed line plots the trend from Arimoto ct al. ( 1996) 
summarizing similar measurements throughout the Local Group. 



recovered in Paper I survey can be accounted for by halos 
of molecular gas around the detected GMCs, comprised 
in part by low mass molecular clouds f M2.3|l . These low 
mass clouds may form in the same event as the high mass 
cloud at the center of the interferometer field. Such low 
mass molecular clouds would have easily been detectable 
in the off fields of the 12-m observation and may be re- 
stricted to regions where high mass molecular clouds also 



form. 

The position angle correlations in i)4.1.4l suggest that 
the formation of molecular clouds occurs on length scales 
less than 500 pc. A common formation event could pro- 
duce aligned molecular clouds, analogous to the forma- 
tion of the solar system producing aligned rotation in 
the planets. Correlations in the orientation of angular 
momentum can be predicted in simulations of turbulent 
cloud formation that include realistic models of galactic 
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shear. Such predictions should show the observed degree 
of alignment with the galaxy and correlations among the 
GMCs. 

The angular momentum problem can be solved by ap- 
pealing to MHD braking. Forming the molecular clouds 
from a narrow range of radii requires a long formation 
time. Formation mechanisms that accumulate gas super- 
sonically cannot have any angular momentum dissipated 
by Alfven waves since the accumulation time is smaller 
than crossing times for magnetosonic waves. In contrast, 
the magnetic field is sufficient to brake forming molecular 
clouds and braking occurs over timescales comparable to 
those of instabilities that could potentially form molecu- 
lar clouds f H4.1.2|l . 

Molecular clouds appear to have their macroscopic 
properties (R, AV) set by their masses (S3- I n addi- 
tion, there appears to be a constant star formation rate 
per u nit molecular mass within galaxies ijWong fc Blitd 
2002). If this is true, then the properties of star forma- 
tion on a galactic scale depend on the amount of atomic 
gas converted into molecular clouds and little else. Al- 
though there are similarities in the properties of molecu- 
lar clouds for a fixed mass, there is good evidence that the 
mass distribution of molecular gas into individual clouds 
varies among galaxies. Different mass spectra may be 
the hallmark of multiple formation mechanisms convert- 
ing atomic into molecular gas. The steep mass spectrum 
for the molecular clouds in M33 (Paper I) represents the 
only difference in the molecular cloud population between 
M33 and the Milky Way. The mass distribution of molec- 
ular gas is thus the dominant factor controlling star for- 
mation on a local and galactic scale. 

5. CONCLUSIONS 

We have presented high resolution follow-up ob- 
servations to a surve y of GMCs in M33 made by 
lEngargiola et alJ ((2003) , using the BIMA array and the 
UASO 12-m millimeter wave telescope. A total of 45 
individual GMCs were detected. The UASO 12-m ob- 
servations measured the total flux in 18 fields, providing 
information on the flux lost by the interferometer and the 
presence of a diffuse component of molecular gas. From 
our observations, we made the following conclusions. 

1 . There is no evidence for a diffuse molecular compo- 
nent traced by CO spanning the disk of M33. We place a 
3(7 surface mass density limit of Sh 2 < 0.3 Mq pc -2 on 
the presence of such a component. This is significantly 
less than the typical surface density of atomic hydrogen, 
7M pc- 2 . 

2. Most of the CO flux in the galaxy is associated with 
GMCs with smaller clouds clustered around the large 
cloud. Low mass clouds may appear in other parts of 
the galaxy but are not ubiquitous. 

3. The velocity gradients of the molecular clouds are 
approximately linear. The magnitudes of the gradients 
are comparable to those found in the Milky Way. If the 
gradients are due to rotation, the rotation period is sig- 
nificantly longer than a cloud lifetime (for t c i ouc i < 30 
Myr). 

4. GMCs show significantly smaller angular momenta 
than are predicted by simple formation theories. The 
discrepancy is, on average, a factor of 5 and ranges up to 



two orders of magnitude. If the velocity gradients are due 
to turbulent motions, the discrepancy widens by at least 
a factor of 2. Both the Toomre and Parker instabilities 
predict angular momentum values that are discrepant 
from observations. 

5. The specific angular momentum is related to mass 
as j(M) cx M°- 6±cu . This value is cons istent with th e 
derived value for the Galaxy of j cx M 0J l|Phillit)slll999|) . 
There is no significant variation of the specific angular 
momentum with galactocentric radius, though such vari- 
ations are predicted by all the large scale accumulation 
theories. 

6. The clouds appear somewhat aligned with the ro- 
tation of the galaxy, though a random distribution of 
position angles is not strongly excluded. If the velocity 
gradients are from rotation, only ~ 60% of the clouds are 
prograde rotators. 

7. The projected velocity gradients of neighboring, 
high-mass clouds are preferentially aligned. This corre- 
lation vanishes for separations larger than 500 pc. This 
alignment may be a signature of the large scale mecha- 
nisms that dictate the formation of molecular clouds and 
should be seen in turbulent cloud formation models. 

8. The observed molecular clouds are similar to those 
found in the Milky Way. We find a size-linewidth rela- 
tionship of AF cx r 0A5:5z0 - 02 and a mass-line width rela- 
tionship of M cx A]/ 4 - 2 ± - 3 . Both of these relationships 
are indistinguishable from those found in the Milky Way. 
There is also no detectable variation in the column den- 
sity of the molecular clouds with the mass or radius of 
the clouds. 

9. Equivalent virial masses and luminous masses 
imply that the X factor for M33 is equal to 2 x 
10 20 H2 cm _2 /(K km s _1 ). There is no significant vari- 
ation of the X factor with metallicity over a range of 0.8 
dex. 

10. These observations support a model for cloud for- 
mation in M33 using atomic gas as the progenitor ma- 
terial for molecular clouds. The formation process accu- 
mulates the atomic gas over a small distance (< 500 pc) 
with significant braking by magnetic fields. 
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APPENDIX 

IMPARTED ANGULAR MOMENTUM 

Consider the progenitor mass for a molecular cloud distributed in galactocentric radius with a linear mass density 
X(R). The specific angular momentum of this material is 



3 = h J x w v w R dR - 



If the cloud's final distance from the center of the galaxy is R Cl then we define r = R — R c and perform a Taylor 
expansion of the velocity curve around R c to get 



3 = jjJ A M [V + vD{R c )\ (r + R c ) dr 



where 



R — R c 



Multiplying out the integral gives, 



1 

3 ~ M 



MVoRc + [R C D(R C ) + Vo] J A(r)r dr + D(R C ) J X(r)r 2 dr 



where the first term represents the angular momentum of a particle orbiting the center of the galaxy and the other 
two terms are the spin angular momentum of the cloud, the measured quantity in these observations. In the absence 
of external forces, the cloud will collapse to the center of mass, which implies / A(r)r = 0. In this case, the imparted 
angular momentum is 

3 gai = V R c + vD(R c )AR 2 

where the matter is accumulated from a region of width 2AR. We define a constant rj set by the moment of the mass 
distribution in terms of the dimensionless parameter u — r/AR: 

r] = J \{u)u 2 du. (Al) 

For collapse in a circular region of radius AR around R c with uniform surface mass density, r] = 1/4. 
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